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Abstract 

Particle physics model building within the context of string theory suggests that 
further copies of the Higgs boson sector may be expected. Concerns regarding tree- 
level flavor changing neutral currents are easiest to allay if little or no couplings of next 
generation Higgs bosons are allowed to Standard Model fcmiions. We detail the resulting 
general Higgs potential and mass spectroscopy in both a Standard Model extension and 
a supersymmetric extension. We present the important experimental constraints from 
meson-meson mixing, loop-induced b ^ sj decays and LEP2 direct production limits. 
We investigate the energy range of valid perturbation theory of these ideas. In the 
supersymmetric context we present a class of examples that marginally aids the fine- 
tuning problem for parameter space where the lightest Higgs boson mass is greater 
than the Standard Model limit of 114 GeV. Finally, we study collider physics signatures 
generic to next generation Higgs bosons, with special emphasis on Ah hhZ 4:b + 2l 
signal events, and describe the capability of discovery at the Large Hadron Collider. 
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1 Generations of Higgs Bosons 

Chiral matter comes in three generations. The simplest hypotheses of electroweak symmetry 
breaking and fermion mass generation assumes the existence of one Higgs boson in the case 
of the Standard Model and a pair of Higgs bosons in the case of supersymmetry. We ask 
here what the consequences are of having more generations of Higgs bosons in analogy to 
fermion matter content. We are not the first to ask this question and investigate answers 
(see for eg. [1], [2] and p] to be highlighted later). Some of our discussion will be known to 
readers, but that is only to set the stage for describing further material we have developed and 
in particular detailing Large Hadron Collider (LHC) implications for next generation Higgs 
boson ideas that survive scrutiny. 

The question is of increased interest of late for two reasons. One, the LHC begins soon 
and enters the prime real estate of Higgs boson phenomenology, and we should be prepared 
to discover all reasonable and viable ideas. The physical particle spectrum of the Higgs 
sector, if it exists, is speculation at present. Investigating various scenarios that may yield 
phenomenology that is different from the simplest Standard Model (SM) approach is needed 
in order to develop more interpretive power over the data when it comes. 

A second reason to consider a next generation of Higgs bosons is from recent developments 
in string phenomenology. Some approaches to particle physics model building from string 
theory suggest that further copies of Higgs bosons may be generic among solutions. For 
example, in the work of [1] a second Higgs generation is generic among the heterotic vacua, 
and may even be more copious than single generation Higgs boson theories. It is also typical 
in this approach that there is a selection rule that allows only the first generation of Higgs 
bosons to couple to the fermions. We will review later why this aspect is very helpful for the 
viability of a next generation of Higgs bosons. 

Next generation Higgs bosons are motivated in other theories as well. For example, in 
theories with branes at singularities bifundamental states come from the same quivers, and 
multiplicities of Higgs pairs are generic just like multiplicities of other representations. In 
intersecting D-brane theories, the chiral content is constrained by topological intersection 
numbers, but the vector-like states can be many-fold. Usually only our self-imposed restric- 
tions in seeking solutions results in one generation. In heterotic orbifold models exotics are 
generic. Restrictions to three families of fermions rarely necessarily restricts Higgs bosons to 
one pair. Some approaches, such as Z3 orbifolds with two Wilson lines [2], naturally provide 
three generations of Higgs bosons, for example. 
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Most physicists nowadays carry the vague suspicions that additional Higgs bosons are 
disastrous unless introduced into very restricted frameworks. They give the photon a mass, 
or result in unacceptable tree-level flavor changing neutral currents. This is to a large degree 
correct, but there are interesting viable limiting cases, touched on above, that are supported 
by theory model building. We set out to elucidate some general conditions for the viability of 
next generation Higgs bosons. We detail a formalism for the analysis, including determining 
the mass matrices and mixing angles in both the SM and Supersymmetry. This culminates 
in a study of a key process at the LHC that is signal for next generation Higgs bosons. 



2 Overcoming Tree-Level Flavor Changing Neutral Currents 

Let us begin by considering an extra Higgs doublet ^extra that is added to the SM Higgs doublet 
$sm- The vacuum expectation values (vev) of each are ($) = Vsm and {^extra) = Vextra, subject 
to the condition that v"^ = v^^ + fg^^^^ = (246 GeV)^. We assume that both Higgs doublets 
couple to the SM fermions. From these two doublets, three degrees of freedom are eaten and 
become longitudinal components of and Zl, and five degrees of freedom are left: the 
scalar mass eigenstates {H,h}, the pseudoscalar A, and the charged Higgs bosons H"^. 

It is always possible to write the Yukawa Lagrangian terms as, 



+V2 CYjQiL^±UjR + V2 CSQiL^±DjR + V2 CfjhL^±EjR + C.C. (1) 

where Ei, Ui and Di are mass eigenstates of leptons, up type and down type quarks and ^^-^ 
are a priori arbitrary. The definition of ^^ev is the linear combination that contains the full 
vev, 

^vev = ^^^^^extra + — *sm, and (2) 
V V 

^_L ^ sm ^extra (3) 

V V 

is the perpendicular state with no vev associated to it. In general, there is nothing to forbid 
the off-diagonal elements of ^fj^'^ from being 0{1). This is the origin of the tree-level Flavor 
Changing Neutral Curent (FCNC) problem of extra Higgs bosons. 

To estimate the experimental upper bound on the off diagonal elements of ^^j^ let us 
assume that the matrices ^'^'^ are real and symmetric. We then obtain the following Feynman 
rules for the scalar and pseudoscalar mass eigenstates (the vertex factor being —i times the 
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Figure 1: Flavor changing neutral current contributions to — mixing from (a) Higgs 
exchange diagrams in an arbitrary 2HDM (there are also t-channel diagrams that we have 
not shown here), and (b) SM gauge contributions. Note that the SM diagrams are one- 
loop whereas the competing Higgs exchange is tree-level. Experiment is consistent with SM 
results, which implies severe constraints on the Higgs flavor-changing neutral current couplings 
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Here the mixing angle a' is the one that rotates from {^yev, to the mass eigenstates {H, h}. 



The most stringent constraints on C,'/j^ come from — mixing (where F = K, B^, D, Bg). 
In a two-Higgs doublet model (2HDM) with arbitrary Yukawa couplings there is a tree-level 
contribution to the F^ — F^ mass splitting because of diagrams like Fig. [l|a). For a' = 
using the expressions derived in we get in the vacuum insertion approximation. 



,2 _|_ 2 



(7) 



where. 



Meson (quarks) 


Bf 


fp (GeV) 


AM|fP' (GeV) 


i^o (ds) 


0.79 ±0.04 ±0.08 


0.159 


(3.476 ± 0.006) X 10"^^ [6] 


(uc) 


0.82 ±0.01 


0.165 


(0.95 ±0.37) X 10-^^ 


SO (db) 


1.28 ±0.05 ±0.09 m 


0.216 ±0.022 


(3.337 ±0.033) x lO'^^ [6] 


B's (sb) 




0.281 ±0.021* 


(117.0 ±0.8) X 10-^3 p 



Table 1: Data associated with the neutral mesons i^", B^ and D^. Values have been obtained 
from [5] unless mentioned otherwise. 

*This is actually the value of the product 



Here Mp is the mass of the meson, niH, rrih and are the masses of H, h and A, Jf is the 
pseudoscalar decay constant and i?^ is the S-parameter of the vacuum insertion approxima- 
tion defined in [1] . We present the values of these parameters and the experimental values for 
AMf in Table [ij 

The Cheng-Sher ansatz |8j is sometimes assumed for the flavor changing couplings, 

^U,D,E X U,D,E ^"^i^i /-.„N 

kij - \j — - — , (lOj 

as many approaches to flavor model building would give rise to it. Let us now find the 
maximum allowed value of consistent with the experimental data and simultaneously 
determine the corresponding ^^/^ ■ For 5° and B^ we use the SM predictions from lattice 
QCD, AM|f = (4.5 ± 1.0) x lO^^^ Q^y |gj ^^^^ AM|f = (135 ± 20) x 10-^=^ GeV [3] and 
add the theoretical error in quadrature to the experimental error. To find the upper bound 
on ^^j^ and C,^ we demand that the sum of the SM value and the 2HDM contribution from eq. 
([T]) does not exceed the experimetal value in Table [l] by more than two standard deviations. 
In obtaining the upper bounds on ^£ and C.^^ we require that just the 2HDM contribution 
from eq. ([T]) does not exceed the experimetal value in Table [l] by more than two standard 
deviations. Taking tuh = mh = mA = 120 GeV we find, 

(11) 
(12) 

(13) 

(14) 

The upper bounds stated above have an uncertainty due to the uncertainty in the value of 
/fVBf- From eq. ([t]) we can see that a 10% uncertainty in fF^/Bp would translate to a 10% 
uncertainty in the upper bound. Greater precision in lattice estimates of the SM values is 
required for more stringent constraints on ^^-^^ from — mixing. 
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Future measurement of branching ratios of rare B decay modes such as B{Bs — fJ'^fi^) at 
LHCb is another way effects of the flavor changing couphngs and may be discovered. 
Expressions for this branching ratio in the SM, and the Higgs contribution in a 2HDM with 
arbitray Yukawa couphngs, can be found in Ref. [10]. The SM value is (3.51±0.50) x 10~^ [10] 
and the current Tevatron upper bound at 95% CL is 5.8 x 10"^ [llj. As LHCb would reach 
the sensitivity to measure the SM value p23, we can estimate the kind of upper bounds 
LHCb experiments would put on and by demanding that the Higgs contribution to this 
braching ratio is less than the SM value. Assuming again that and are real and equal 
to each other, and that ^ij = Xij^/mirrij/v, we obtain 

A£ < 2 ^^"^^ ^ GeV) (Tevatron inferred bound) (15) 

A£ < 0.6 ^""-^/^^Q ^^^^ (LHCb expected sensitivity). (16) 

Although this limit from Tevatron and future sensitivity expectations at LHCb do not appear 
to be as powerful as meson-meson mixing constraints, the uncertainty in what A^^ might be 
in the presence of next generation Higgs bosons suggests that it should stay under consid- 
eration. For a more thorough discussion of B^ — B^ mixing constraints and the processes 
Bs — !■ — Kjj,^ji~ see Ref. [TU]. For a discussion on the constraints on the diagonal 

couplings 4^ see Ref. [13] . 

We have seen above that the off-diagonal couplings are required to be extremely small 
in order to satisfy FCNC constraints. There is a general class of solutions to this problem [T3] 
while admitting the existence of extra Higgs bosons in the spectrum. Tree-level FCNCs do 
not arise if Higgs boson interactions with the fermions take the form 

^Cf = ygQ[^ F,({$,}) D;.^ + F„({<|.,}) f/j^ + yf<L[^ F,i{^,}) E'^^ + c.c. (17) 

where all the quark fields are gauge eigenstates. Fu^d,e{{^k}) are functions of Higgs fields 
constrained only by the requirements that they are independent of the fermionic flavor 
indices i,j and that transforms like an SU{2)l doublet with hypercharge —1/2, and Fd 
and Fg transform like SU{2)l doublets with hypercharge 1/2. 



The generalized form of eq. ( 17) subsumes many ideas already present in the literature. For 
example, the SM Higgs sector is F^ = H^m and Fd = F^ = Hsm- The type II [I5] 2HDM [I6] 
is Fu = Hu and Fd = Fe = Hd- The type I 2HDM [16] is F„ = ^\ and Fd = F^ = $i with 



an additional $2 that does not couple to fermions. The leptophilic Higgs model of Ref. [TTj is 

Fu = Fd = (pq and F^ = (f)i. 
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There are an infinite variety of models tliat can satisfy eq. (17). However, principles are 
expected to be at work to fall into this class if there is more than one Higgs boson. In the 
case of supersymmetry, the type II structure follows from holomorphy of the superpotential. 
In the case of type I theories, it is usually assumed that the second Higgs has, for example, a 
discrete Z2 symmetry ^2 = —^2 that forbids its direct coupling to fermions whereas $1 does 
not. 

The summary point is that a next generation Higgs boson is unlikely to satisfy FCNC 
currents due to tree-level mediated interactions unless a principle is invoked the ensures that 



it will satisfy the condition of eq. (17). The most straightforward principle that we can invoke, 
and one that has been nicely illustrated recently in the model of [T], is a selection rule that 
forbids the next generation Higgs boson from coupling to fermions. 

3 Next Generation Higgs bosons of Supersymmetry 

We begin with a discussion of next generation Higgs bosons in supersymmetry. In minimal 
supersymmetry there are already two Higgs doublets present in the spectrum. In unrestricted 
field theory, two Higgs doublets with arbitrary couplings are a disaster for flavor changing 
neutral currents. However, as mentioned above, supersymmetric theories have the special 
property that all superpotential operators must be holomorphic in the superfields. Thus, it is 
impossible to employ one Higgs field to give masses to both up-type fermions and down-type 
fermions. The introduction of the second Higgs doublet solves this problem, and holomorphy 
is the principle by which tree-level FCNC's are held under control. The interaction lagrangian 



takes the form of the type II 2HDM in the form of eq. (17). Of course, when supersymmetry 
is broken, non-holomorphic interactions can induce additional couplings leading ultimately to 
potentially interesting flavor changing neutral currents [12] , but those are naturally small over 
much of parameter space. 

If we wish to add more Higgs doublets to supersymmetry, we must do it in up- and 
down-Higgs boson pairs. This is in sympathy with adding a new generation - a new copy - 
of the B^dS pair. More importantly, it is required to straightforwardly satisfy anomaly 
constraints due to the presence of new fermions that are the superpartners of the Higgs bosons. 
When a next generation of Higgs doublets is added to the spectrum, holomorphy is no longer 
powerful enough to save us from tree-level FCNC. Additional Yukawa couplings generally 
create FCNC challenges [21 [13]. In Ref. [3] supersymmetric three generation Higgs models 
have been considered where an exact symmetry prevents the next generation Higgs Bosons 
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from coupling to the fermions. They show that these next generation Higgs bosons would 
not couple to the standard Higgs bosons and under the assumption of equal scalar masses 
at unification scale they would not get vevs, thus avoiding FCNC problems. In general the 
simplest way out of the FCNC challenge is to assume that the extra Higgs doublet pairs couple 
very weakly, or not at all, to the fermions. This is the assumption we shall adopt for now. 

The result of the discussion above is that we are overlaying a type I Higgs structure to our 
type H supersymmetric theory. In other words, we are adding a Higgs boson pair that does 
not couple to the fermions on top of a Higgs pair that does. Our emphasis in this study is on 
the type I aspect (i.e., the next generation), and as such the supersymmetric type II feature is 
of less immediacy. Thus, we shall postpone a detailed discussion of additional constraints to 
later, to the less complex model of adding one Higgs boson to the SM. Nevertheless, we wish 
to present the mass eigenstates and mixing, and in particular highlight what effect multiple 
Higgs bosons have on the mass of the lightest Higgs boson scalar of supersymmetry. 



3.1 General Higgs Potential 

For the supersymmetric two Higgs generation model (SUSY 2HGM), we consider two addi- 
tional Higgs superfields Hu2 and Hd2 charged under SU{2)l x f/(l)y just as Hui,Hdi respec- 
tively. The terms in the superpotential involving these superfields are, 

W — HiiHulHdl + fJ,l2HulHd2 + IJ'2lHu2Hdl + //22-f^u2-f^d2- (18) 

From now on by Hui, Hdi, Hu2, Hd2 we will denote only the scalar part of the superfields which 
have the components. 

The Z^-terms are given by, 

= -g Y,{Hi;r^H^i + HlfT^Hd^- (20) 

i 

The D-term contribution to the potential is given by, 

a 

i i 
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where A 
bution. 



-. Using the form of the superpotential in eq. (18) we get the F-term contri 



Vf = (l/^llP + \^ll2?)\Hul\^ + (|^22r + \^l2l?)\Hu2? + (l/illT + \^2l?)\Hdl\^ 

+ (|/X22r + \lJii2?)\Hd2? + {auHl^Hu2 + adHlHd2 + c.c.) (22) 



where, 



O-d — /^ll/^12 + /^2l/^22 
O'u = /^ll/^21 + /^l2/^22- 



(23) 
(24) 



In addition to the D-terms and F-terms there are the soft SUSY breaking terms, 



hi2HuiHd2 + h2iHu2Hdi + h22Hu2Hd2 + c.c.) 



where H,,jH, 



Adding all the terms from eq. (21), (22) and (25) we get finally, 



V 



(25) 



Vd + Vf + Vsoft 

ml^\Hui\'^ + ml^\Hu2? + ml^\Hdi\'^ + ml^\Hd2\^ + {auHl^Hu2 + adH\-^Hd2 + c.c] 

+ {bilHulHdl + hi2HulHd2 + h2lHu2Hdl + fe22-f^«2-f^d2 + c.c.) + \^[^{\Hui[ 



\Hdn' 
(26) 



Here is the sum of and the F-term contribution from eq. (22). All the couplings are 



required to be real by hermiticity with the exception of the bilinear couplings hij and Oj which 
are in general complex. By redefining the phases of the doublets Hd2 and Hu2 only three 
of these six couplings can be chosen to be real so that the theory is in general CP violating. 

Note that we are allowing Hui^i to mix arbitrarily with Hu2,d2 in the Higgs potential only. 
Some mixing between the two generations is necessary as we do not want to have a global 
symmetry under which the next generation Higgs bosons can be rotated with respect to the 
first generation because this would lead to the presence of a massless Goldstone boson. When 
it comes to interacting with the fermions, as discussed before and emphasized again later, next 
generation Higgs bosons are generally barred from having couplings unlike the first generation 
Higgs bosons. 
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3.2 Electroweak symmetry breaking and scalar mass matrices 



We put (-ff^) = (-f^rfj) = 0, {H^) = -^Vi and demand that the first derivatives of the potential 
with respect to the fields Hf vanish to obtain, 



m 



^n^ - h."^ + + A(E-i - E<) = (27) 



dl 



Vdl Vdl Vdl 

J J 



- i.„^ - bn^ + -aj^ - v% - J: = (28) 

•^■ul '-'ul 'Jul 

3 J 



m 



6..^ - + + HJ:-1 - E<) = (29) 



d2 



m 



Vd2 Vd2 Vd2 

Vd2 7 Vdl ~ Vul ,,>r^ 2 2 



u2 



h,^ - k,^ + -au-^ - X{Y.-l - E<) = 0- (30) 



Vu2 Vu2 Vu2 

3 3 



Here hij = Re bij and Oj = Re Oj. It is necessary that (Hl^j) = {H^^ = for electromagnetism 
to remain unbroken. We can always choose one of the charged fields, say H^^ to have an 
expectation value (-f^^^i) — using the SU{2)l gauge freedom. To ensure that {H^i) = 
= (i?^) = is consistent with the minimization conditions, however, we must demand 
that in addition to eqs. (27)-(30) the second derivatives of the potential with respect to the 
charged fields at this point are positive. This is equivalent to demanding that the masses of 
the three physical charged Higgs bosons are positive. 

Also note that we have assumed the vev of the neutral components Vui,Vdi to be real and 
positive. One of the vevs say, Vdi, can be chosen to be real and positive using the f/(l)y 
gauge freedom. As we discussed earlier this theory is in general CP violating therefore we can 
choose the other vevs v^i, Vu2 and Vdi to be real and positive simply by a convenient choice 
of phases of the respective doublet fields. These phases can then be absorbed in the bilinear 
couplings bij and and a redefinition of the quark fields. If the underlying theory ensures 
that all the bilinear couplings bij and are real so that the Lagrangian conserves CP, in order 
to avoid spontaneous CP violation Vui, Vu2 and Vdi must be real too. Therefore the point 
{Aui) = {Ad2) = {Au2) = must be a minima where Ai = Im (i/f). For this, in addition 
to eqs. (27)-(30) the second derivatives of the potential with respect to the pseudoscalar fields 
Ai must be positive also at this point. This is equivalent to demanding that the three physical 
pseudoscalars in the theory must have positive masses. Once Vu2 and Vd\ are known to be 
real they can always be chosen to be positive by a convenient choice of signs of the respective 
doublet fields, which can then be absorbed in the the bilinear couplings bij and Oj and a 
redefinition of the quark fields. If bij and aj are complex the CP even states mix with the 
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pseudoscalar states. To avoid this unnecessary complication to our present purposes we will 
present below the mass matrices assuming that the bij and are real. 

Using eqs. (27)-(30) we can eliminate the mf. We can find the mass matrix of the four 
CP-even scalars and the four pseudoscalars by substituting Hf = + hi + iAi). The mass 
matrix for the CP-even scalars, M'jj, in the {di,ui,d2,U2} basis is, 



( 

-bn - 2\vuiVdi 
ad + 2XvdiVd2 

\ -621 - 2\Vu2Vdl 



-fell - 2\vuiVdi 

2 

-bi2 - 2\VulVd2 
+ 2XVulVu2 



ad + 2\vdiVd2 

-bi2 - 2\VulVd2 



m 



H33 



-621 - 2\Vu2Vdl \ 

au + 2At;„it;„2 

-622 - 2Xv^2Vd2 



J22 



2\Vu2Vd2 



m 



I 



where. 



U22 



'f^Hu = (^n- ^^21;;^ ad- h2At;rfi 

+ 2At;^i 

b^9— + bm— - ad— + 2Af, 
621 — + bm— - a,„ — + 2\v. 
The pseudoscalar mass matrix, A^^, in the {^1,^1,^2,^2} basis is. 



m 



H33 
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H4A 





Vul 


+ O21 

Vdl 


Vd2 


11 


Vdl 


- ad — 

Vdl 
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+ O12 

Vul 


Vu2 


11 


Vul 
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Vul 


, ^"2 
+ 622 — 
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Vdl 


12 


Vd2 


- ad — 

Vd2 




Vdl 


+ O22 

Vu2 


Vul 


21 


Vu2 


au 

Vu2 



2 

d2 



2 

u2- 
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2 
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where, 
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A22 
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Vd2 

I — 

Vdl 
Vu2 



m 



A33 



bu h 621 — 

Vdl Vdl 
, Vdl , , Vd2 

bn h O12 — 

Vul Vul Vul 

, Vul , , Vu2 Vdl 

bi2 h 622 ad — 

Vd2 Vd2 Vd2 

, Vdl , Vd2 Vul 
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The charged Higgs mass matrix differs from the pseudoscalar mass matrix only due to the 

+ V:,o - v:jo Vd^ Vu^ Vd^ Vdl f rfl Vu2 \ 



last term in the potential in eq. (26). In the basis {H^* , H^^, H^^ 1 ^^2} ^ 



Ml 



Ml 



+ 



3_ 
4 



V 



Ml + ^«2 ~ '^d2 
VulVdl 
Vd2Vdl 
Vu2Vdl 



'dl 



VdlVul 
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3.3 Upper bound on the mass of the hghtest CP even Higgs 



We want to transform the mass matrices above to the Runge basis which is defined as follows. 
One of the basis vectors in the Runge basis is, 

Vi = Vdl/V + Vul/V H^i + Vd2/V H^2 + Vu2/V Hy_2- (31) 

Here v — -x/f^i + f^i + + '^d2- ^di (2) l conjugate representation and is given 

by, 

Hi^(^). (32) 

We choose the other basis vectors so that they are orthogonal to this vector and to each other. 
The simplest choices for two of the basis vectors are, 

V2 = Vy,i/vi Hl^ - Vdi/vi Hui 

— * 

Vs = Vu2/V2 H^2-^d2/v2 Hu2 (33) 

where Vi = a/^^^]~+^^ ^2 = \/v^2 + '^d2- ^® fourth basis vector which is 

orthogonal to the first three by the expression, 

VU = {V^.U)V^ + {V2.U)V2 + {V3-U)V3 - U 

V'a 

V, = , ' (34) 



where U can be any arbitrary vector. The transformation matrix is an 5'0(4) rotation matrix, 

R=iy^V2V3Vi). (35) 

In the Runge basis, the mass matrices of the pseudoscalar, charged Higgs and CP-even 
scalar Higgs boson masses take on special form. For the pseudoscalar and charged Higgs mass 
matrix, the first row and first column contain all zeros, which is expected since the first basis 
vector Vi is the "electroweak vev multiplet" which has all the vev. Thus, the CP-odd and 
charged components of the first basis vector in the Runge basis are the spin-zero Goldstone 
boson states absorbed by the W^, and W£ vector bosons. What remains is one more 
component of a full doublet, namely the CP-even part. We shall write the full CP-even mass 
matrix below and identify the matrix element corresponding to the mass of this CP even scalar 
and comment on its meaning. 

The CP even mass matrix in the Runge basis is, 

= R^MlR. (36) 
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0.0 0.2 0.4 0.6 0.8 1.0 

sin 00 



Figure 2: Contours of the upper bound on the mass of the hghtest Higgs (in Mz units) in the 
SUSY 2HGM taking cos 2/32 = — 1- We have also shown contours of equal A^^ (in GeV), the 
energy scale at which the top Yukawa coupling At becomes larger than Att. Details about the 
calculation of Ksc appear in section [5j 
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where, 



m 



SXiVulVdl)'^ + buvf + b2lvl-^Vu2 - duV^iVul - advliVd2 + buVdiVd2 



'2 

H22 — 2 



ViVulVdl 



'2 _ adVulVu2 + bl2VdlVu2 + b2lVulVd2 + auVdlVd2 + 8XVulVu2Vd2 
^ ViV2 

'2 _ v'^{b2iVuiVu2 + auVdiVu2 " adVuiVd2 - bi2VdiVd2) + 8X{yuiVdivi2 - vl^Vdivl^) 

mj^24 — 2 

VV{V2 

,2 _ 8X{Vu2Vd2f + b22V2 + &211'dlWrf2 ~ auVulVd2 " + ^12WwlfL 

"^//33 — 2 

V2Vu2Vd2 

'2 _ v'^{b2lVdlVd2 - auVulVd2 + adVdlVu2 - bl2VulVu2) + 8X{vj^vl2Vd2 - vliVu2Vd2) 
f^HM — 2 



VV1V2 

2 „,2 „,2 „,2 ^2 



V {b2lVdlVu2 - auVulVu2 + bi2VulVd2 " adVdlVd2) + 8A(V^^<2 - ^ul'"d2) 



^Hii ~ 9 2 2 



The Runge basis helps us see what the hghtest Higgs boson mass becomes in the hmit that 
supersymmetry breaking masses are large, fh ^ Mz- In that case, the {11} element of the 
CP even Higgs boson mass matrix is the only diagonal element that stays small. A theorem 
of linear algebra tells us that the smallest eigenvalue of a positive definite matrix is smaller 
than the smallest diagonal element. Therefore, at tree-level we know from this {11} element 
the upper bound on the mass of the lightest CP even Higgs h, which i^ 

< ^2 

=^ rrih < M^l cos^ cos 2/3i + sin^ cos 2/32I, (37) 
where tan/3j = Vui/vdi and tanw = V2/V1. 

The leading supersymmetry breaking corrections to this expression are from top squark 
loops in the same manner as found in the Minimal Supersymmetric Standard Model (MSSM). 
Thus, the controlling difference between our theory and the MSSM prediction for the Higgs 



mass is the tree-level expression of eq. (37) compared to that of the MSSM, which is 



mh < Mz\ cos2P\. (38) 

One sees that if we set tan /3i = tan /3, which becomes in both theories the fixed value for the 
ratio of vevs of the two Higgs doublets that couple to the fermions, a small additional contri- 

^Notc that such an upper bound would exist even if our assumption that hij and are real is not true. 
For complex bij and the CP even and pseudoscalar states mix. There are again seven neutral scalars and a 

/ 2 ^2 \\2 

Goldstone boson. The 7x7 mass matrix of the physical scalars would again have ^"^^-^"1?-""'^^ as the {11} 
element in the seven dimensional basis in which Vi is the Runge vector. Here the Runge vector in the original 
eight dimensional space is Vi — 'i;2, "^3, 114, 0, 0, 0, 0). 
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sinu mt (TeV) Asc (GeV) 



2.7 2 X 10" 

0.5 2.0 5 X 10* 

0.7 1.4 3 X 10^ 

0.9 0.8 2 X 10^ 



Table 2: For various values of sinw we show the value of rht needed to raise the Higgs boson 
mass above the experimental limit of 114 GeV and also the scale Asc where the top Yukawa 
coupling diverges. The tree-level contribution is obtained by assuming tan /3i = 1.5 and 
cos 2/32 = -1- 

bution can be made to the Higgs boson mass in our next generation Higgs theory compared 
to the MSSM by allowing for larger tan/32 > tan/3i = tan/3. 

The available gain to the Higgs boson mass in this manner is tiny if tan /3i > 5. For 
lower values of tan Pi the additional Higgs doublet pair contributions can be significant if the 
mixing angle u is larger. This may be useful since it is a challenge in the MSSM to obtain a 
Higgs boson mass above the 114 GeV experimental limit without having too-high superpartner 
masses that induce fine-tuning in the electroweak sector potential. 

In Fig. [2] we plot contours of the tree-level Higgs boson mass in units of Mz in the plane 
of tan/3i vs. sincu, assuming that cos 2/32 = — 1? which is a good approximation if tan/32 ^ 5. 
One sees that as sin co increases, the second Higgs doublet is more responsible for electroweak 
symmetry breaking and therefore the Higgs mass increases due to larger tan ^2- The drawback 
is that the Higgs bosons that couple to the fermions get smaller vevs, leading to larger Yukawa 
couplings. The larger top Yukawa coupling could diverge at a low scale. We describe this 
effect in more detail in section |5j In Fig. |2] we have also plotted, therefore, the contours of 
the scale Asc at which the top Yukawa coupling becomes strongly interacting (i.e., diverges). 
In Table |2] we show that with increasing sinu mixing angle, it is possible to have a smaller 
value of the stop masses so that that the Higgs boson mass is greater than 114 GeV. The 
quoted values of rht = y/rht^rht^ are obtained by assuming that the tree-level contributions 
are derived from tan/3i = 1.5 and cos 2/32 = —1, and that only the leading order radiative 
correction is contributing to the Higgs boson mass. 



For higher values of sino;, the tree-level contribution increases, thereby putting less pressure 
on the stop masses to raise the Higgs boson mass above 114 GeV. This is clear from the values 
in the table, where for higher sino; the needed rht values are lower. 




(39) 
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We also show in Table [2] the scale Asc at which the top Yukawa coupling diverges for the 
various values of sincj. As is expected, the larger the value of sinw for some given tan/3i, 
the lower the Asc scale. This is the tradeoff one has between a larger Higgs mass prediction 
and a lower scale of top quark Yukawa coupling divergence. This is reminiscent of the Next- 
to-Minimal Supersymmetric Standard model (NMSSM) which can have an arbitrarily large 
Higgs boson mass by adjusting the parameter A in the superpotential interaction XNHuHd, 
however, large values of A imply divergences well below a putative unification scale. This is a 
common feature in many attempts to solve the Higgs mass bound challenge of supersymmetry. 



4 Next Generation Higgs boson of Standard Model 

4.1 Electroweak symmetry breaking and scalar mass matrices 

To add a next generation Higgs boson to the SM is equivalent to postulating a 2HDM with 
two scalar doublets $i and $2 having hypercharge 1/2. Earlier we discussed the many ways 
that a second Higgs boson can be added to the spectrum without being incompatible with 
experiment. There are many options, including type I models and type II models and variants 
on that theme. The type II structure is most naturally incorporated within supersymmetry 
using holomorphy as the guiding principle, as we discussed in the previous section. Going 
beyond that, the most straightforward way to incorporate extra Higgs doublets is to implement 
a type I structure. In other words, the second Higgs boson (or next full generation) induces 
no tree-level FCNC by virtue of it having no Yukawa couplings with the SM fermions. 

To ensure no couplings of the second Higgs boson to fermions, the discrete symmetry, 

$2 ^ -$2. (40) 

can be imposed, for example. If we allow a soft violation of this symmetry by dimension-two 
terms we can still avoid tree level FCNC bounds. The most general renormalizable potential 
for the scalars in which the discrete symmetry is softly broken only by dimension-two terms 



V{^i,(^2) = /i|$2r + &($l$2 + c.c.) + Ai|<l>ir + A2|$ 

+ A3|<l'in<l'2|^ + A4(<I>5$l)($I$2) ■ 



^($t$2)^ + C.C. 



(41) 



The potential above is the potential in the toy model of [Ij once the singlet in their theory 
gets a vev. Hermiticity requires all the coupling constants in the above potential to be real 
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with the exception of b and A5. A convenient choice of the phase of $2 will allow only one of 
the two couplings to be real so that the theory is CP violating in general. Note that if both b 
and A5 vanish the potential is invariant under a global U{1) symmetry for $2- Thus either b 
or A5 must be non-zero to prevent the pseudoscalar from being a massless Goldstone boson. 

Two CP conserving limits of this potential have been considered in the literature [19]. 
The limit in which 6 = has been called potential A (Va) and the limit in which A5 = has 
been called potential B (V^). Let us now analyze in some detail the electroweak symmetry 
breaking (EWSB) pattern in these two limit^ 

Potential 6 = 0, A5 7^ 

Potential A can be obtained by strictly imposing the discrete symmetry ^2 ^ —^2 which 



requires 6 = in eq. (42) 



K4(<i'i,$2) = + /i'i<i'2r + Aii$ii^ + A2i$2r + A3i$iri$ 

+ A4($5$l)($I$2) + 



^($l$2r + C.C. 



(42) 



Without loss of generality all the Aj couplings are real. Hermiticity demands it for all Aj 
except A5, which can always be rotated to real and chosen to be either positive or negative by 
making $2 absorb its phase. For definiteness we choose A5 < here. 

The potential must be bounded from below in all field directions. One can test for danger- 
ous runaway directions by parameterizing field excursions arbitrarily large in value. The fol- 
lowing field directions give us all the unbounded from below (UFB) constraints (see, e.g., [21J): 

$f , ^2 direction UFB constraint 

(0,a), (0,0) Ai > 

(0,0),(0,a) A2>0 (43) 

(0,A^/'a),(A;/V0) A3 + 2v/A^>0 

(0, Xl^\), (0, X^^a) A3 + A4 + A5 + 2^/XJr2 > 0. 

The most general vacuum expectations values for the two $1^2 Higgs fields can be expressed 
as (see, e.g., [2^ ) 

A non-zero U2 would indicate the full breaking of SU{2)l x U{1)y, and in particular the 
photon would obtain mass. For electromagnetism to remain unbroken when $1 and $2 get 



^For a discussion about the vacuum structure and the possibihty of spontaneous violation of CP and 
electromagnetism in more general 2HDMs see Ref. [20] , 
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vevs the following condition must hold (see for instance |22l E]), 

A4 + As < 0. (45) 
The minimization conditions obtained by setting dV/d(j)i = for all real fields 0j defined in 

$i=f t'^'t'], and $,= f t'^'t'], (46) 
then ensure that EWSB is proper and the doublets get the vevs, 

*' - ( JV2 ) • *^ - ( Jv-2 ) <^^) 

where f 1 > and we can always choose f 2 > by a convenient choice of sign of the doublet 
field $2- The minimization conditions dV/dcf)^ = and dV/dcj)^ = given by |22l 

can be used to eliminate the parameters fii and fi2 |23j . 

We need to check if this solution is stable. To do that we require that the second derivative 
of the potential, i.e. the mass matrix, be positive definite. We can find the mass eigenvalues 
are by solving four 2x2 matrices. These matrices arise from (f)k4'k+A mixing for = 1, 2, 3, 4. 



We simplify the entries in these matrices by substituting nf and from eq. (48). To begin 



with, we look at the 0i05 and 0206 mixings, which have the same 2x2 mass matrix: 

which leads to four eigenstates 

m^± = (charged Goldstone bosons) (50) 
m^± = — 'i^i+'^l) (charged Higgs bosons). (51) 

The mixing angle is. 



tanw = — . (52) 



Now let us look at 0408 mixing: 



-^*4<^8 - 1 x,v^V2 -X^vl ' • ^^^^ 



19 




Figure 3: Parameter space for massive and massless photon in a type I 2HDM with potential 
A. 

This leads to two eigenstates 

rriQ = (neutral Goldstone bosons) (54) 
m\ = —X5{vl + V2) (neutral pseudoscalar boson). (55) 

The mixing angle is again tan uj = V2/V1. Finally, there is 0307 mixing: 

..2 _ ( 2Ait;? (A3 + A4 + \-,)viV2 \ 

-^<^3<^r - ^ (A3 + A4 + \,)vrV2 2\2vl J ■ 

This is the 2x2 mass matrix for the two physical neutral scalar Higgs bosons of the theory, 
h and H. The mixing angle to rotate from {03,07} basis to {H, h} basis is usually called a, 
which is defined by convention to satisfy 



H \ f cos a sin a 
h j \ — sin a cos a 



(57) 



The solutions are obtained by simple eigenvalue, eigenvector analysis of the 2x2 matrix, and 
one obtains 



tan 2a = - — 2 ^ — 2 (^^ 



A345^lf2 

where A345 = A3 + A4 + A5. We will use the convention that H is the 'SM-like' Higgs and hence 
we will always take a < 7r/4. The masses of the CP-even scalars are. 



^hH = Ai^^r + ^2vl ± V (Aif? - \2vlY + X^ii^vivl. (59) 
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We require that the mass matrix be positive definite, which puts important constraints on 
the parameters of the theory. For example, from the charged Higgs and pseudo-scalar Higgs 
boson masses we know that 

A4 + As < 0, and A5 < (60) 



is required. Note that the first condition is the same as eq. (45), the condition that ensures 
that the photon remains massless. Fig. |3]plots the parameter space in the A4 vs. A5 plane that 
corresponds to massive photon and massless photon cases in agreement with Ref. [22] • 

Potential B: A5 = 0, 6^0 

Potential B can be obtained by imposing a global U{1) symmetry for $2 and allowing it 
to be broken only softly by dimension- two terms like 6$J$2; 

1/($1, $2) = + /i'|$2|' + h{^\^2 + C.C.) + Ai|$i|^ + A2|$2|' + A3|<I>in$2r 

+A4($^$i)($l$2). (61) 

Without loss of generality all the Aj couplings are real due to hermiticity. The coupling h can 
be rotated to real by $2 absorbing its phase. 

The conditions for the potential to be bounded from below in all field directions are the 



same as eq. (43) with A5 = as these conditions are determined by the quartic couplings so 



that the additional bilinear term does not affect them. The most general vacuum expectations 



values for the two $12 Higgs fields can again be expressed by eq. (44) 



The minimization condition dV/d(j)i = [22] in this case is given by, 

bu2 + ^U2f if 2 cos ^ = (62) 

It is clear that for 

& + -^fif2 cos^ 7^ (63) 



it is required that U2 = from the 0i minimization condition in eq. (62). This ensures that 
electromagnetism is not broken and the photon remains massless as the doublets get vevs. 
This condition and the other minimization conditions then ensure that the doublets get vevs 
of the form, 

*' = ( vjv2 ) ■ ^""^ *^ = ( vJv-2 ) 
where f 1 > and we can always choose f 2 > by a convenient choice of sign of the doublet 
field $2; which can then be absorbed in a redefinition of the coupling b. The minimization 
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conditions dV/dip^ = and dV/d(f)i = given by |22j . 



03 : 1^1 + bv2 + ^^vl + Xivl = 
07 : + bvi + ^^vf + X2VI = 

can be used to eliminate the parameters /xi and fi2 |23j . 



Let us now look at the mass matrix after substituting fif and from eq. (65). Let 
look at the 0i05 and 0206 mixings, which have the same 2x2 mass matrix: 

which leads to four eigenstates 

mQ± = (charged Goldstone bosons) 
v? + v^ / 2b 



mfr± = I A4 H (charged Higgs bosons), 

V1V2' 



the mixing angle being. 



, V2 

tan uj = — . 



bv2/vi b 
b -bvi/v2 



Now let us look at 0408 mixing: 

This leads to two eigenstates 

~ (neutral Goldstone boson) 

m\ = ("^1 + "^2) (neutral pseudoscalar boson), 

V1V2 

where the mixing angle is again given by tanw = V2/V1. Finally, there is 0307 mixing 

2\iV^-bV2/Vi {\s + \i)ViV2 + b 

{Xs + ^4.)viV2 + b 2X2vl-bvi/v2 



P7 



We have again. 



H \ ( cos a sin a 



^3 



h j \ — sin a cos a. j \ (l)^ 
The masses and mixing angle are, 

A34f 1-^2 + b 



tan 2a 



XxVx - MV2 + 
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Figure 4: Parameter space for the different patterns of EWSB in the type I 2HDM with 
potential B. 

where A34 = A3 + A4, and 



m 



I ^ = \^vl + \2vl - ± J (\,vl - \,vl + fe^l-^) ' + {\,,v,v, + hf. (76) 

2^ViV2 V V 2,ViV2 J 



Requiring the charged and pseudoscalar Higgs mass to be positive gives us, 

A4 + — < — < 0. (77) 

V1V2 V1V2 

Fig. |4] shows the the regions in the parameter space which lead to the different patterns of 
EWSB for potential B. 



An analysis of the general potential in eq. (41) is complicated by the fact that both h and 



A5 cannot be chosen to be real by redefining the phase of $2- This leads to a mixing between 
the pseudoscalar and CP even states. If both h and A5 are assumed to be real, however, the 



masses and mixing angles for the general potential in eq. (41), assuming that the vacuum is 
invariant under CP and electromagnetism, are, 

2 ^? +^2 A ^ 26 \ 

= A4 + As + 78 

m\ = -{vl + vl)(\, + —\ (79) 



<,H = ^ivl + ^2vl - 6^1^ ± y [Xivl - \2vl + fe^l;^) ' + A345t^l^i (80) 

tanw = — (81) 



Vl 
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Figure 5: The solid lines show contours of constant r, the fractional deviation of the type I 
2HDM value of the 6 — )■ 57 width from the SM value. The dashed lines are the boundaries 
of the region in the m^f-f-tana; plane excluded by constraints from the b s'j branching 



rate at two standard deviations. As one can see from eq. (84) there are two different ways of 



satisfying the constraint: (1) if tanu is small the type I 2HDM width is close to the SM value, 
or (2) tan^ w can be tuned to a higher value so that J2i=ct ^iCnimb) in eq. (84) is close to 



the SM value {Cn ~ — (7^*^). Thus there are two disconnected allowed regions in Fig. [s} The 
star shows our choice of and tanw that we use for collider simulations later. 



tan 2a 



Xivf - \2vi + b^ 



(82) 



where A345 = A3 + A4 + A5. Thus in 2HDMs from the original eight degrees of freedom we 
get three Goldstone bosons and five physical scalars. The Goldstone bosons are absorbed as 
longitudinal modes by W^l~ and Z. The tree level masses of the W~^^~ and Z are, 



Mw = Mz cose 



w 



^3) 



We now give some of the constraints that this model experiences when requiring compat- 
ibility with all past experiment. 



4.2 Indirect Constraints 

Constraints due to virtual H^/^ effects 

In the type I 2HDM the coupling of H^/~ to the fermions is proportional to tanw. Thus 
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in the type I 2HDM, unlike the type II 2HDM, all the constraints coming from processes 
involving virtual if"*"/" can be met for small enough tanw. The strongest such constraint 
comes from the fe — )■ 37 branching rate. An expression for the width of the process in type I 
2HDMs is [21], 

2 



Cn{Mw) 



1287r4 

l = C,t 

V-''/''[Cn{M^) - - 1) - ^(r^^^/^^ - 1)] 



Here Aj = V*.Vib, t] = as{mh) /as{Mw), Aw is the SM contribution due to W^^^ loops and 
Ah is the additional contribution due to triangle diagrams involving the charged Higgs. For 
expressions of Aw and Ah see for instance [21]. The SM value for the width can be obtained 
by simply putting u; = or taking the limit mH+ — )■ 00 in the above equation. 

Let the fractional deviation of the type I 2HDM value from the SM value be, 



^2HDM-I 
b—^s-f 



-pSM 



^5) 



Fig. [5] shows the contours of constant r in the m/^+-tanco' plane. The dashed lines in Fig. |5] 
show the region excluded by experiments at 2a level. To obtain the dashed curves we have 
used the world average for the experimental value of the branching fraction [25j . 

B{B X,7) = (3.55 ± 0.24t°:°^ ± 0.03) x 10"^ 



We have rescaled the theoretical value obtained using eq. (84) to reproduce the NNLO SM 
prediction [26] . 

B{B X,7) = (3.15 ± 0.23) x 10"^ 
in the m/^+ — )■ 00 limit and added the error associated with this value in quadrature to the 



experimental error. From eq. (84) one can see that there are two different ways of satisfying 
'1) if tanw is small the type I 2HDM width is close to the SM value, or (2) 



the constraint: 
tan^ u can be tuned to a higher value so that 



m eq. 



84) is close to its SM 



value {C-ii ~ — Cfj*^). Thus there are two disconnected allowed regions in Fig. [s} We can see 
from Fig. |5] that the constraint is satisfied for all values of mH+ if, 



which translates to, 



tanw < 0.32 



V2 < 75 GeV. 



?6) 
^7) 
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For 'mH+ = 200 GeV we obtain the constraint tanw < 0.47 {v2 < 105 GeV). Another con- 
straint due to virtual H^^~ effects comes from — oscillations. This, however, puts a 
weaker constraint than the b ^ sj process at the 2a level 



It is important to note that a small V2 not only implies a small u but also suggests a 



small a for a 2HDM with the general potential in eq. (41). This can be seen from the 



expression for tan 2a in eq. (82). Using eq. (79) we can obtain the expression for tan 2a for 
tanw = V2/V1 ^ 1, 

tan 2a ^ ^ , , J^' ^\ tan uj 88 

which shows that a small f 2 suggests that a should be small too. 

The p 'parameter 

In the 2HDM the p-parameter gets additional contributions from corrections to Mpv and 
Mz due to loops of the scalars. With respect to the SM theory with a Higgs boson mass of 
rrih = 120 GeV, the range of Ap that can be tolerated p] by replacing the single SM Higgs 
boson with the 2HDM is 0.0000 < Ap < 0.0012 at the 68% CL, where Ap = p{2HDM) - 
pijnfj^^). The 2IIDM computation needed for this can be found in Ref. [28]. One finds that it is 
easy to satisfy precision electroweak constraints for a 2HDM with masses in the neighborhood 
of 50 — 300 GeV, as we shall see later in our example benchmark points that we use to study 
collider signatures. 



4.3 Collider Constraints 

The fact that none of the scalars h,H or A were seen at LEP puts constraints on their masses 
mh, mn, t^a and 5 = |a — We are especially interested in the cases when the decay modes 
A — i- Zh 01 H ^ hh/AA are kinematically allowed, i.e. when the conditions — frihl > mz 
or rriH > 2mh/A hold respectively. 

The major constraint on the mass of the non SM-like CP-even Higgs, rrih comes from the 
non-discovery of h produced by the Higgstrahlung process e^e~ — )■ hZ at LEP. The cross 
section of the Higgstrahlung process in the 2HDM is suppressed by a factor of sin^ 6 with 
respect to the cross section of the same process in the SM. It is therefore possible to have h 
as light as we wish by choosing a sufficiently small value of 6. Upper bounds on the cross 
section for this process for a particular value of mh therefore give upper bounds on sin^ 6. 
Fig. |6] shows the LEP2 upper bounds [29j on sin^ 6 vs m^ assuming h decays entirely to bb. 
As one can see from the figure there are no constraints on mh at all, if sin 5 < 0.1, and no 
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Figure 6: Upper bounds on sin^ 5 as a function of nih derived from searches for the Hig- 
gsstrahlung process e~^e^ — > Zh at LEP with center of mass energy ^/s = 91-209 GeV. It has 
been assumed that h decays entirely to bb. We also plot the three parameter sets A, B and C 
that we have chosen for simulations (Table |3| on the sin^ 6 - rrih plane. The curve has been 
reproduced from Ref. [29] . 

constraints on sin 5 if nih > 114 GeV. If we take sin 5 < 0.2, nih ~ 100 GeV is safely within 
allowed limits. 

Another important process that could have been potentially seen at LEP is the associated 
production process e+e^ — hA . The cross section for this process is proportional to cos^ 6 
in the 2IIDM. Analysis of LEP2 data, assuming — mh\ > so that A predominantly 
decays as A — t- Zh, has been done by the DELPHI collaboration [30] and the process has been 
found to be unconstrained. Even if this condition does not hold there are no constraints on S 
if we have + m,, > 200 GeV 

The constraints on the pseudoscalar mass ttia are much weaker. The cross section of 
associated production process, e"'"e^ — )■ HA, is proportional to sin^ 6. The results of LEP 
analyses motivated by this process put upper limits on sin^ 6 for a given and mfj (see 
Ref. [29] for details). If the dominant decay mode is — ?■ AA there are no constraints at all 
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if niH > 120 GeV and tua > 50 GeV [29]. 

As far as the charged Higgs is concerned, the LEP2 direct search constraints from the 
process e~^e~ — ?■ Z* — )■ H^H~^ place a lower limit of 76.7 GeV on [21]. The Tevatron 
search for t — t- H^h puts a 95% CL upper bound on B{t — > H^h) at 0.1 — 0.3 in the mass range 
90 GeV< mH+ < 150 GeV assuming B(i7+ — ?■ cs) to be 100% [32j. For almost the entire mass 
range B{t — )■ H^h) < 0.1 is allowed. In the type I 2HDM this branching ratio is proportional 
to tan^ uj, so that we obtain tanw (~ s^) < 0.3 for mH+ = 90 GeV and tan w < 0.8 (s^ < 0.6) 
for mH+ = 150 GeV by requiring B{t — t- H^b) < 0.1. The limit on tana; is actually too 
conservative for m/^+ = 90 GeV as the upper bound on the branching fraction for this mass 
is about 0.3. 

Finally let us consider the SM-like Higgs boson H. If itlh < 2m;i/^, H would predominantly 
decay into quarks with branching ratios very similar to that in SM. In this case the process 
e'^e'^ —7- HZ puts a weaker lower bound on mn in 2HDMs than the SM value 114 GeV because 
the 2HDM cross section of this process is smaller than the SM cross section by a factor cos^ 6. 
li H hh/AA is allowed there are of course no constraints on 6 if ttlh > 115 GeV, and for 
6 > 0, rriH can be even smaller [29]. We will take mn = 120 GeV for all the parameter sets 
we use in our simulations in section [6l 

5 Yukawa Coupling Perturbativity 

If there exists at least one Higgs boson with a vev that does not couple to a fermion /, 
the Yukawa coupling of that fermion Xf must necessarily be greater than its corresponding 
would-be SM value. Larger Yukawa couplings in the theory run the risk of renormalizing 
to strong coupling at a lower scale than desired. In our study this is a consideration that 
must be explored, since our emphasis is on next generation Higgs bosons that do not couple 
to the fermions. Therefore, in this section we quantify where Yukawa couplings blow up in 
renormalization group evolution as a function of the vev of fermiophobic Higgs doublets in 
the theory. The effect is particularly pronounced for the top quark Yukawa coupling, since it 
is associated with the highest mass fermion in the theory. 

The one loop Renormalization Group Equations (RGBs) for the Yukawa (A/) [33] and 
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Figure 7: Asc, the energy scale at which either one of the three large Yukawa couplings A^, A;, 
or A^ becomes larger than An, as a function of c^. We show curves for the type I 2HDM 
case as well as the the SUSY 2HGM with different values of tan/3i. For the tan /3i = 50 
and tan /3i = 65 cases we show two dashed curves corresponding to the minimum (left) and 
maximum (right) allowed value of mb(M^) in eq. (102). For the mean value mb{Mz) = 2.85 
GeV in eq. ( 102[ ), for tan/3i > 55 it is the bottom Yukawa which becomes strong (> Att). We 
have used the Renormalization Group Equations at one loop level. 



gauge couplings (see for eg. [MJ) in an extension of the SM with rid extra Higgs doublets is, 
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+ 3 YliD -^D + YliE ^ ^^'^ ^ denote the up-type quarks, the 

type quarks and the leptons while g' , g2 and g^ are the U{1)y, SU{2)l and SU{?>) 
couplings. Note that the summation over U, D does not include a summation over 



down- 
gauge 
colors. 

There are similar equations for the other up- type quarks (t — )■ f/ in eq. (89)), the other down- 
type I 



type quarks (6 — )■ D in eq. (90) and the other leptons (r — )■ in eq. (91)). For the 
2HDM Ud = 1. 
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For supersymmetric models with Ug extra generations of Higgs doublet pairs than MSSM 
the one loop RGBs for Yukawa [33] and gauge couplings (see for eg. [35]) are, 
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For the supersymmetric two Higgs generation Model (SUSY 2HGM) of section 3.1 
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There are similar equations for the other up-type quarks (t — )■ f/ in eq. (93)), the other 



down-type quarks (6 — )■ -D in eq. (94)) and the other leptons (r — )■ in eq. (95)) 



The Yukawa couplings in type I 2HDM are given by, 

Xfiruf) = V2- 
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where = Vi/v (with v = 2Miv/g). In the SUSY 2HGM the Yukawa couplings are given by, 

(98) 
(99) 
(100) 

where tan/3i = Vui/vdi- 

In the type I 2HDM all the Yukawa couplings except for the top Yukawa coupling can be 
ignored. In the SUSY 2HGM in addition to the top Yukawa, the tau and bottom Yukawa 
couplings also become important at high values of We have solved the RGE with the 
boundary conditions, 

Aasirrit)' 
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2.7-3.0 GeV [36] 
1.75 GeV 
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Set Input Parameters Ap 



A 




= 120 GeV, ruh = 


= 50 GeV, 


rriA 


= 150 GeV, 


0.0010 




mH+ 


= 200 GeV,s^ = 


0.1, Sa = 


0.2 1 


smo = O.lOj 




B 


rriH 


= 120 GeV, ruh = 


= 70 GeV, 




= 180 GeV, 


0.0003 




mH+ 


= 200 GeV,s<^ = 


0.1, Sa = 


0.2 ( 


sin (5 = 0.10) 




C 


rriH 


= 120 GeV, rrih = 


100 GeV 




= 200 GeV, 


< 0.0001 






= 200 GeV,s^ = 


0.1, Sa = 


0.3 ( 


sin 5 = 0.20) 





Table 3: Example parameter sets. The values of Ap are computed in these 2HDMs with 
respect to the SM value with Higgs mass of 120 GeV. 



and as{Mz) = 0.1182. In eq. (102) the allowed range corresponds to the Particle Data Group 
(PDG) range for mb{mb) = 4.03 — 4.37 GeV p|. Let Asc be the strong coupling scale, i.e 
the energy scale at which either one of the three large Yukawa couplings A^, Af, or At becomes 
larger than Arc. In Fig. [7] we show how Asc varies as a function of c^^. We show the curve for 
the type I 2HDM case as well as the curves for the SUSY 2HGM for various values of tan 
In the SUSY 2HGM at high values of tan /3i ^ 50 the bottom Yukawa coupling becomes strong 
at a lower energy scale than the top Yukawa coupling. 

We see from the results of this section that if we allow the fermiophobic next generation 
Higgs doublet too large of a vev, the enhanced Yukawa couplings that are required to make 
up for the smaller vev of the Higgs boson that the fermions couple to may diverge at a lower 
scale than desired. For example, if one wishes to preserved supersymmetric gauge coupling 
unification up to the scale of ~ 10^^ GeV there are critical values of cos uj that cannot be 
crossed depending on tan/3i, which leads to maximum values of the next generation Higgs 
boson vev. 



6 Signatures at the Large Hadron Collider 

In this section we identify processes that can provide signatures of next generation Higgs 
bosons at the LHC. While our numerical results have been obtained for the type I 2HDM our 
basic conclusions are true for next generation Higgs bosons in general. 
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6.1 Dominant decay modes 



We will now compute the branching ratios of the various decay modes of the neutral Higgs 
bosons in the type I 2HDM. The relevant Feynman rules can be found in Ref. [T9] . 



Fig. [8[a) shows the branching ratios for decay of the pseudoscalar A (see eq. (42)). The 
decay modes A — )■ hh, A — )■ HH, A — )■ WW and A — )■ ZZ are not allowed by symmetry as 
A is a pseudoscalar. As shown in the figure, when allowed kinematically, the branching ratio 
for A — )■ Zh is nearly unity. 

The SM-like Higgs H decays mainly via the modes H — )■ hh/AA even if modes like H — )■ 
WW and H ZZ are kinematically allowed. Fig. [sj^b) shows the branching ratios of H in 
potential A with a light h. 

The branching ratios of h are very similar to that of a SM Higgs for small values of 
when decay modes like h — t- ZA and h — )■ A A are not kinematically allowed. For rrih > 2171 a, 
the decay mode h — )■ AA overwhelms all other modes including h — )■ ZA, h — )■ WW and 
h — )■ ZZ. A very interesting limit is a — )• 0, w — 0. In this limit h becomes both fermiophobic 
and bosophobic (the tree level coupling of h to fermions is proportional to —Sa/cuj and the 
coupling to vector bosons is proportional to sin 5) and the dominant decay mode becomes 
h — )■ 77. This case has been dealt with in detail in Ref. [12] and [TF] . 

6.2 The pp — > Zh{bb)h{bb) signal and choice of input parameters 



As we noted in section |4.2 indirect constraints put an upper bound on V2 in the type I 2HDM 



and this not only implies a small co but also suggests a small a and hence suggests a small 
6 = \a — u\. A small vev for the fermiophobic Higgs doublets is also required for high scale 
perturbativity as we saw in the last section. If 6 is small a light h can satisfy the constraints 
due to the e+e^ — )■ hZ at LEP (Fig. |6]). The pseudoscalar A must then be chosen sufficiently 



heavy to satisfy the LEP constraints from e^e hA (see section 4.3). 



Another implication of a small 6 is that unlike the process pp hZ, which is highly 
suppressed (by sin^5 in 2HDMs), the processes pp — )■ Ah and pp — ZH are only mildly 
suppressed (by cos^ 6 in 2HDMs). In the first case, if allowed kinematically, A predominantly 
decays as A — t- Zh (Fig. [sj^a)) and in the second case H predominantly decays as if — )■ hh/AA 
(Fig. ^h)) if allowed to kinematically so that both processes can lead to the same final state 
Zh{bb)h{bb)F 



Note, however, that m a 2HDM with potential B, in the hmit a,iuj — > we get — ?> mj^ (see eq. (72) 



and eq. (76l) so that A — > Zh would not be kinematically allowed in this limit. This issue does not arise for 
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Figure 8: (a) Branching ratios of the various decay modes of the pseudoscalar A with mu = 120 
GeV, rrih = 70 GeV, mH+ = 150 GeV, Sa = 0.2 and s^^; = 0.1 in any 2HDM and (b) branching 
ratios of the SM-hke Higgs H in potential A with = 180 GeV, rrih = 70 GeV, mH+ = 150 



GeV, 



0.2 and s.., = 0.1. 



Table [3] shows three different sets of input parameters for which we perform simulations to 
compute the pp — )■ Zh{bb)h{bb) signal. Set A receives contributions from both the pp — )■ Ah 
and pp — )■ ZH processes whereas all the other sets receive contributions from the pp — )■ Ah 
process only. We have also computed Ap for all the data sets in Table |3] to show that this 
constraint is met. As for the other constraints, we have marked these parameter sets in Fig. |5] 
and Fig. |6]to show that our parameter sets satisfy the 6 — 37 constraint and the e^e~ — > hZ 
constraint respectively. 



As mentioned in section 4^ a light A is experimentally less constrained than a light h. 
Although we will not perform simulations for the process pp — )■ A{bb)h followed hj h ^ ZA{bb) 
the analysis of the process would be very similar to the process we will consider. This process 
would be important if /i — t- ZA is kinematically allowed but h AA is not kinematically 
allowed. If mn > 2mA the process pp — ?■ ZH{AA) — ?■ Zbbbb will be a very important 
signature. Although we will not perform a simulation for such a scenario, we will compute the 
contribution to the cross section of the process pp — )■ ZH{hh) — )■ Zbbbb for Set A (this will be 
a part of the net cross section). This is again expected be very similar to the case when A, 
rather than h, is the lighter scalar that H decays into. 

It should be clear that all the arguments that have led us to the pp A{Zh)h 



potential A or the general potential in eq. ( 41 1 . 
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Figure 9: The 2b invariant mass spectrum Mf,i, for the signal on top of the SM background 
for the input parameters in Set A (see Table [s]), showing a peak at ruh = 50 GeV. Two 
processes contribute to the signal (1) pp — > A{Zh)h — )■ Z{l^l^)Zbbbb and (2) pp — > ZH{hh) — )■ 
Z{l^l^)Zbbbb . The cuts applied are those in eqs. (104)-(107) and the center of mass energy 
has been taken to be 14 TeV. The branching ratio B{H — )■ hh) ^ 1. The reconstruction 
efficiency of the lepton pair and the jets and the 6-tagging efficiency have all been taken to be 
unity at this stage. 



Zh{bb)h{bb) signal are valid not only for the type I 2HDM but for any multi-Higgs theory 
having scalar doublets that do not couple to fermions. It is easiest to satisfy the constraints 
from charged Higgs loop contributions (as in the 6 — )■ 37 process) and high scale perturbativity 
if the vevs of these fermiophobic doublets are small. Let h and A be mass eigenstates that 
contain mostly the CP-even and CP-odd neutral components of such a doublet respectively. 
A small vev of this doublet would imply that the ZZh coupling strength is small, but the 
ZAh coupling strength would not get suppressed. 



6.3 Signal and Background cross section at LHC 

We now present the results of the simulations we performed for the pp — i- Zh{bb)h{bb) signal. 
The analyses we present are new; however, there are related studies that we will point out 
to the reader in our Set A analysis. All our results are independent of the choice of 2HDM 
potential except for the H ^ hh contribution to the signal in parameter set A that depends 
on B{H — hh) which is model dependent. B{H hh) has been taken to be equal to its 
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Figure 10: The 46 invariant mass spectrum M^f, for the signal on top of the SM background 
for the input parameters in Set A (see Table [3]). Two processes contribute to the signal (1) 
pp A{Zh)h Z{in-)Zhhhh and (2) pp ZHQih) Z{l+l~)Zbbbb . The pp ZHQih) 
mode gives rise to a peak at mn = 120 GeV. We have smeared the 46 invariant mass assuming 
an experimental resolution of 20% of mn (24 GeV in this case) for the reconstructed peak. 
The cuts applied are those in eqs. ([I04l)-([l07l) and the center of mass energy has been taken to 



be 14 TeV. The branching ratio B[H — hh) ^ 1. The reconstruction efficiency of the lepton 
pair and the jets and the 6-tagging efficiency have all been taken to be unity at this stage. 

value in potential A which is nearly 1. We used MADGRAPH [38] to generate signal and 
background events at 14 TeV center of mass energy for the process pp Zh{bb)h{bb) for 
the different parameter sets in Table |3} We then decayed the Z to l^l^ (/ = e, /i) using the 
DECAY software in the MADGRAPH package. Note that we have ignored any contribution 
to the SM background from final states with lepton pairs not produced in Z decay that have 
invariant mass close to Mz nevertheless. The following basic selection cuts have been applied 
using MadAnalysis, 

(104) 
(105) 
(106) 

where AR = (A?])^ + (A^)^ and p is any parton (i.e. quark or lepton) in the process. 
The background cross section with these cuts is 9.9 fb. To reduce the background further we 
use the fact that two 66 pairs have the same invariant mass in the signal up to experimental 
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Figure 11: Invariant mass spectrum M^b for the signal on top of the SM background for the 
input parameters in Set B (see Table |3]), showing a peak at rrih = 70 GeV. The signal process 
is pp — 7- A{Zh)h — 7- Z{l^l^)Zbbbb . The cuts applied are those in eqs. (104)-(107) and the 
center of mass energy has been taken to be 14 TeV. The reconstruction efficiency of the lepton 
pair and the jets and the 6-tagging efficiency have all been taken to be unity at this stage. 



resolution. Detector simulations of scalars decaying into 6-pairs [32] find that because of 
detector effects like calorimeter energy resolution, electronic noise, and physics effects like 
final state radiation, energy loss outside cone and semi leptonic decays, only about 85% of 
the events register di-6-jet invariant masses within 20% of the true value. To simulate this 
effect we smear the invariant masses according to a Gaussian distribution such that 85% of 
the events lie within 20% of the mean. There are three ways to divide the four 6-quarks 
(say abed) into two pairs {ab cd, ac bd and ad be). The combination that gives the invariant 
masses of the two pairs (after smearing) closest to each other has been considered (note that 
in reality the experimental uncertainties in the invariant masses of the six possible 6-pairs are 
not uncorrelated as assumed here). Let M^bi and Mbb2 be the smeared invariant masses of 
the two b pairs thus selected, and let Mbb be the mean of these two numbers. We impose the 



following cut in addition to those in eq. ( 104 ) 



\Mbbi-Mbb2\< 0.2 Mbb. (107) 



36 



■ Signal 




60 80 100 120 140 160 180 

Mbb(GeV) 



Figure 12: Invariant mass spectrum M^i, for the signal on top of the SM background for 
the input parameters in Set C (see Table |3]), showing a peak at nih = 100 GeV. The signal 
process is pp — )■ A{Zh)h — i- Z{l^l^)Zhhhh. The cuts applied are those in eqs. (104)-(107) and 
the center of mass energy has been taken to be 14 TeV. The reconstruction efficiency of the 
lepton pair and the jets and the 6-tagging efficiency have all been taken to be unity at this 
stage. 



With this additional cut the background is reduced from 9.9 fb to 2.6 fb whereas the signal 



is only reduced to about 70% of the value with only the cuts in eqs. (104)-(106). We provide 



the cross sections for the different parameter sets and the background in Table |4j Note that 
the background cross section can be further reduced by requiring to be in a certain mass 
window around uih 

In Fig. [9] we plot the 26 invariant mass spectrum for Set A obtained for the events 



passing the cuts in eqs. (104)-(107). For Set A the contribution to the cross section of 



pp — > Zh{bb)h{bb) comes from two different processes pp — )■ A{Zh)h — > Zbbbb and pp — 
ZH{hh) — Zbbbb. The branching ratio B{H — hh) has been taken to be equal to its poten- 
tial A value which is nearly 1. As shown in Fig. [9] the contribution from pp — )■ ZH(hh) mode 
is quite large. This contribution also gives rise to a peak in the 46 invariant mass spectrum as 
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Set 



Cross Section 
with the cuts 



Cross Section 
with the cuts 



in eqs. (104)-(106) (fb) in eqs. (104)-(107) (fb) 



Detector level Cross- 
Section (cTe// in eq. (108)) 
with the cuts 



in eqs. (104)-(107) 



Set A 

pp — >■ Ah 
pp — 7- ZH 
Total 

Set B 

pp Ah 

Set C 

pp —J- Ah 

SM Background 



1.2 
3.6 
4.8 



0.9 

0.6 
9.9 



0.8 
2.7 
3.5 



0.6 

0.4 
2.6 



0.1 
0.4 
0.5 



0.09 

0.06 
0.4 



Table 4: Signal and background cross sections for pp — )■ Zh{hh)h{bb) at LHC. The center of 
mass energy has been taken to be 14 TeV and the acceptance cuts are those mentioned in the 
top row. 



shown in Fig. [10| We have smeared the 46 invariant mass assuming an experimental resolution 
of 20% of tuh (24 GeV in this case) for the reconstructed peak. We expect very similar results 
if A rather than h is the scalar H decays to. The signal from pp — )■ ZH{AA) Zbbbb will 
be a very promising signature for the type I 2HDM, especially because the constraints on a 



light A are rather weak as discussed in section 4.3 Similar analyses for H — )■ AA have been 



done before (in Ref. [10] the final state lA{bb)A{bb) has been considered while in Ref. [H] the 
final state W{EJp^^'^l)A{bb)A{bb) has been considered). These papers, however, do not apply 
the cut in eq. ( 107[ ) which leads to improved significance of the signal. 



Fig. 11, [12] shows the 2b invariant mass spectrum for the other parameter sets. Only the 
process pp — )■ A{Zh)h — )■ Zbbbb contributes to the cross section in these cases. 

To get the cross section we expect the detectors to effectively measure we must multiply 
by the efficiency of reconstruction of a lepton pair and that of four jets. These efficiencies 
depend on kinematical quantities like pt and 77. We take an average value 0.8 for lepton pair 
reconstruction efficiency (see pgs. 72-92, pgs. 210-223 in Ref. [12]) and 0.9 for reconstruction 
efficiency of a jet (see pgs. 286-287 in Ref. [12]). We also require that at least three of the jets 
are 6-tagged which gives an overall 6-tagging efficiency equal to {^el{l — e^) + where the 
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6-tagging efficiency for single jet tb ~ 0.5 |12]. Putting it all together we get, 

a^ff = 0.16 a. (108) 

This equation is applicable for both the cross section and the background. The effective signal 
and background cross sections also appear in Table |4} 

We have not computed the contribution due to mistagging of c-jets or other light jets. This 
will decrease the signal significance somewhat as mistagging is expected to have an appreciable 
contribution only to the SM background and not the signal. This is because the signal cross 
section is proportional to the square of the branching ratio of h to the mistagged quarks but 
the branching ratio of h to quarks other than 6-quarks is much smaller. In Ref. [H] the signal 
and background cross sections for the process pp — )■ WH — )■ W{EJp'^^'^l)A{bb)A{bb) have been 
computed including mistagging effects. A rough estimate of the cross sections they obtain can 
be made from Fig. 5 in their paper. Such an estimate shows that the background cross section 
they obtain due to mistagged quarks is about a third of the contribution due to correctly 
tagged b quarks. This mistag background is highly dependent upon the details of detector 
performance issues that will be sorted out in the course of the LHC runs. We do not expect 
the addition of these backgrounds to substantively change the discovery capability that we 
have presented, especially since they are unlikely to peak at the Higgs mass rrih in the di-jet 
invariant mass spectrum. 

7 Conclusion and Discussion 

In conclusion, we have argued that next generation Higgs bosons should be viewed as generic 
possibilities in string theory model building, and illustrated this viewpoint with recent devel- 
opments in string phenomenology. We presented a generalized theorem for the structure of 
Higgs couplings to SM fermions that automatically avoids problematic tree-level flavor chang- 
ing neutral currents that are induced by new Higgs boson exchanges. Our viewpoint is that 
the interaction rules of this theorem are too restrictive to be satisfied without a principle. 
In the case of the two-Higgs doublet model of supersymmetry, the principle is holomorphy. 
Additional Higgs bosons added in any other context, such as more Higgs doublet pairs in 
supersymmetry or simply another Higgs boson in the SM, requires a strong discrete symme- 
try or selection rule. This can be contemplated within effective field theories, for example 
by $ — 7- — $ ^2 symmetry, or within string theory by algebraic selection rules that may not 
totally forbid the unwanted couplings but can approximate zero, as was the case in the work 

aim- 
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Next wc considered various constraints that these theories must face. For example, al- 
though tree-level flavor changing neutral currents induced by neutral Higgs exchanges may 
be satisfied, loop-level ones induced by the charged Higgs boson may not be. The 6 — )> 57 
transition is quite constraining to exotic Higgses that get large vacuum expectation values, 
since they steal vev from the Higgs that couples to fermions, thereby raising those fermion 
Yukawa coupling magnitudes to dangerous levels. Even normal Cabibbo-Kobayashi-Maskawa 
(CKM) mixings can create too-large amplitude shift in that case compared to what experi- 
ment allows, and therefore the parameter space is not completely open and limits are derived 
on the exotic Higgs vev as a function of the charged Higgs boson mass. 

We computed the spectrum of Higgs boson states with a next generation in the super- 
symmetric and non-supersymmetric context. Within supersymmetry we showed that there is 
the prospect of slightly raising the tree-level CP even Higgs boson mass with respect to the 
MSSM in the smaller tan/3 region. The effects are largest when the next generation Higgs 
boson has a large vacuum expectation value. The large Yukawa couplings that are present 
when an exotic, fermiophobic Higgs doublet takes a large vev can alter the domain of pertur- 
bativity of the theory. We showed both in the SM context and the supersymmetry context 
that the top quark Yukawa coupling could develop a Landau pole well below the Planck scale. 
A low-scale Landau pole would preclude the existence of a perturbative theory description of 
unification at the high scale, and these results must be taken into account when considering 
a next generation of Higgs bosons. 

Finally, we investigated the phenomenology of the exotic Higgs sectors at the LHC. Multi- 
Higgs boson phenomenology within supersymmetry is a mature topic; however, the fermio- 
phobic next generation Higgs boson layer to the theory has not been considered in depth. The 
salient new features are similar to SM phenomenology with an additional fermiophobic Higgs 
doublet. Thus, we discussed collider physics possibilities within that less complex framework. 

A particularly interesting possibihty is the production of hA at the LHC, followed by 
A — > hZ. Assuming h is rather light, say less than ~ 150 GeV, we can expect the largest 
branching fraction of h decays to be to bb. For good distinction from background we can 
also demand the Z decay to leptons Thus, the signal becomes Ab + 21. Background 

becomes particularly limited when we require at least three 6-quark tags and that the four 
jets reconstruct two equal mass resonances. We propose that as a search strategy for this case, 
and show that there are good prospects for the LHC to find this signal. Discovery would be 
an indication of next generation Higgs bosons. 
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